Macroscopic quantum tunneling in multigap superconducting Josephson junctions: 
Enhancement of escape rate via quantum fluctuations of Josephson-Leggett mode 
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We theoretically study the macroscopic quantum tunneling (MQT) in a hetero Josephson junction 
formed by a conventional single-gap superconductor and a multigap superconductor such as and iron- 
based superconductors and MgB 2 . In such a Josephson junction more than one phase difference 
is defined. We clarify their phase dynamics and construct a theory for the MQT in the multigap 
Josephson junctions. The dynamics of the phase differences are strongly affected by the Josephson- 
Leggett mode, which is the out-of-phase oscillation mode of the phase differences. The escape 
rate is calculated in terms of the effective action renormalized by the Josephson-Leggett mode at 
zero-temperature limit. We successfully predict drastic enhancement of the escape rate when the 
frequency of the Josephson-Leggett mode is less than the Josephson-plasma frequency. 
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Macroscopic quantum tunneling (MQT) 14 is a coun- 
terintuitive phenomenon in quantum mechanics appear- 
ing in a macroscopic level and has been observed in var- 
ious fields of physics such as condensed matters, nuclei, 
cosmology, etc. This phenomenon has still attracted a 
great interest in physics communities. In particular, the 
MQT in Josephson junctions, which is observed in a 
switching event at low temperature— i 4 -, has been inten- 
sively studied because it is promising for applications to 
a Josephson phase qubit 5 7 . 

In this paper we explore the physics of MQT in an un- 
explored type of Josephson junctions which has multiple 
tunneling channels. Such a Josephson junction can be 
fabricated by using recently discovered iron-based super- 
conductors—"— or MgB 2 11 13 , because these supercon- 
ductors are multiband ones having more than one dis- 
connected Fermi surfaces and the superconducting gap 
can be individually well defined on each Fermi surface. 
In a Josephson junction made of multigap superconduc- 
tors one may expect that the superconducting tunneling 
current has multiple channels between the two supercon- 
ducting electro desi 4 ^!. We construct a theory for the 
quantum switching (i.e., MQT) in Josephson junctions 
with multiple tunneling channels. To our knowledge no 
theory has been formulated for the MQT in multigap sys- 
tems. The theory predicts that the escape rate, i.e., the 
rate of quantum tunneling, is drastically enhanced com- 
pared with that in conventional single-channel systems. 

In multigap superconductors a collective mode called 
the Leggett's mode — ~— appears in the low energy re- 
gion, which is an out-of-phase oscillation mode of the su- 
perconducting phases. In Ref. 14 a theory for the Joseph- 
son effect in superconducting hetero junctions made of a 
single-gap superconductor and a two-gap one is formu- 
lated. In such Josephson junctions, because two kinds 
of gauge- invariant phase differences can be defined, there 



are two phase oscillation modes, i.e., the in-phase mode 
and the out-of-phase one, which correspond, respectively, 
to the Josephson-plasma and the Josephson-Leggett (JL) 
mode. In this paper we construct a theory for the MQT 
in superconducting hetero junctions, incorporating the 
degree of freedom of the JL mode into the quantum 
switching event from no n- voltage to voltage states. It is 
shown that the zero-point motion of the JL mode signifi- 
cantly enhances the MQT escape rate when its frequency 
is less than the Josephson-plasma frequency. We also 
point out that the ratio Ej/E[ n in addition to Ej/Eq 
governs the boundary between the classical and quan- 
tum regimes, where Ec, Ej and E- in are, respectively, the 
charging energy, the Josephson coupling energy between 
the two superconductors and the inter-band Josephson 
coupling energy in the two-gap superconductor. 

Consider a hetero Josephson junction made of a single- 
gap superconductor and a two-gap one 14 i 16 i 17 ? as shown 
schematically in Fig.[TJ Such a junction has been already 
fabricated, using the multigap superconductors MgB 2 21 
or iron-based superconductors 22 1 23 . In this system one 
can define two gauge- invariant phase differences, 0^ and 
0( 2 \ Then, the Josephson current density between the 
two superconducting electrodes is given by the sum of the 
superconducting currents in the two tunneling channels 
as ji sin^ 1 ) + 32 sin#( 2 ), where ji is the Josephson criti- 
cal current density in the ith tunneling channel. When a 
voltage v appears between the two superconducting elec- 
trodes, the gauge- invariant phase differences show tem- 
poral evolution satisfying the generalized Josephson re- 
lation 14 
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with oti = e/ii/d and A = 1 + 0102/(01 + 02), where e 
is the dielectric constant of the insulator with a thick- 
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ness d and is the charge screening length due to the 
electrons in the ith band. The constant is related to 
the charge compressibility in the two-gap superconduct- 
ing electrode—. 

As shown in Ref. 14 , the Lagrangian in the hetero 
Josephson junction with a in-plane area W and capac- 
itance C = eW/Aitd is expressed as 

L= o 7^2 ~a l~ ■ - V' + Ej— 0,(2a) 
2 (2e) 2 y A a 1 J ra 2 J l c 

V = -E n cos (9 (1) - £ J2 cos (9 (2) - ^ in cos 0, (2b) 

under a bias current J e x in the absence of an external 
magnetic field, where and are the center-of-mass 
phase difference and the relative phase difference defined 
as 

= _22_0(i) + _^_0( 2 >, v = ^ (1) -^ (2) - 

Ql + CK2 Oil + Q^2 

The first two terms in Eq. (|2b|) are the Josephson cou- 
pling energies with the coefficients Eji = hWji/2e and 
the third term represents the inter-band coupling en- 
ergy, where the coefficient E- in is expressed as E- in = 
HW\ Jin|/2e 14 . Since the "inter-band current" J\ n can 
take both signs, depending on the gap symmetry, we in- 
troduce the sign factor n = J- m /\J' m \. The total critical 
Josephson current I c and the coefficient Ej in the last 
term in Eq. (|2aj) are defined as I c = W\ji + ft^l and 
Ej = HI c /2e. We note that the voltage v appearing in 
the junction is related to only 0, as seen in Eq. (pQ). 

From Eq. (|2a|) one can derive the Euler-Lagrange equa- 
tion for the center-of-mass phase difference 6 as 

A-^ + c^sinflW+a&sin** 2 ) = (3) 

with fi^ Pi = ^/2EcEji and fi^jp = \/2EqEj. We note 
that cjpi is the Josephson-plasma frequency in the ith 
tunneling channel. From Eq. (|2a|) we also have the Euler- 
Lagrange equation for i/j as 

$ + ^cjj L sin = -OLiUJ 2 n sin + a 2 cjj 2 sin # (2) , (4) 

where cjjl is the frequency of the JL mode 14 given as 
Tw^jl = + <^2)^c^in- The above two equations 

are coupled since 6^ and 0^ are functions of 6 and -0. 
We note that the bias current is the source for the time 
evolution of 6 but not for 0, which is consistent with 
the generalized Josephson relation (pQ). It should be also 
noted that we have two characteristic energy scales, the 
Josephson-plasma frequency ujji and the JL one cjjl, in 
this system. 

Let us now study the macroscopic quantum effects in 
the Josephson junction with multiple tunneling channels 
on the basis of the Lagrangian (|2a|) and evaluate the 
MQT escape rate. In the following we assume n > 0, 
since the case of k < shows qualitatively no difference. 

Suppose that the switching to the voltage state is in- 
duced by the quantum tunneling of the phase differences 
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FIG. 1: (Color online) Schematic view of a SIS hetero Joseph- 
son junction. We have two tunneling channels between the 
two superconductors with the critical current densities ji and 
j2 as indicated in the right panel. In the upper two-gap elec- 
trode the inter-band Josephson coupling with the coupling 
constant Ji n exists. 



0^ and #( 2 ) which are confined inside a potential well. 
When both 9^ and 0^ show the tunneling at the switch- 
ing, its transition probability is given by the expectation 
value of the time evolution operator with respect to the 
state |0W = 0, 0^ = 0)(= |<9 = 0, = 0)) 3 , which yields 
the formula for the MQT escape rate as 

r = ^Im^({0},{0};/3). (5) 

Here, the symbol {0} means (9,ip) = (0,0), and j3 is 
the inverse temperature, j3 = l/k&T. The propagator 
K{X, X'\ j3) in Eq. ([5j) is expressed in terms of the imag- 
inary time path-integral 

r X(h/3)=X 

K{X, X , ;f3)= VOV^j e~ fo P d - L */\ 

JX(0)=X' 

where X = (0, i/j) and L E is the Euclidean version of the 
Lagrangian (|2a|) . Let us assume that ip is confined in a 
small region around = at the tunneling, which is 
justified when the inter-band coupling is not so strong. 
In this case one can utilize the expansion with respect to 
-0. Then, up to the order of 2 the Euclidean Lagrangian 
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FIG. 2: (Color online) Schematic energy diagram for the fully 
"quantized" system with two quantum variables and 0. In 
the case where is weakly oscillating within a potential well, 
the energy levels of coincide with those of a harmonic oscil- 
lator with frequency cjjl- The energy levels of is corrected 
by the quantum oscillations of 0. 
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FIG. 3: (Color online) Ratio r/r in the case of Ej /E c = 10 2 . 
We assume that ai = ot2 — 0.1 and ji = j2 for simplicity. 
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FIG. 4: (Color online) Renormalized potential V cm ,ef£- The 
potential for is modified by the zero-point fluctuation of the 
JL mode. 



L E is approximated as L E 
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4(ai + a 2 )-E c \<W '2 
^tat =g+Ejip 2 cos# - g-Ejtpsm6. 



E in i> 2 , (6b) 
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Here, Ej = Eji + Ej 2 and A « 1 is assumed. The cou- 
pling constants g+ and g- in Eq. (f6c|) are defined as g + = 
(^ji/2^j)[ai/(ai+a 2 )] 2 + (^ J2 /2Ej)[a 2 /(ai+a 2 )] 2 and 
0- = (Ej 1 /E J )[a 1 /(a 1 + a 2 )} - (E J2 /E J )[a 2 /(a 1 + a 2 )]. 
We note that in the fully quantum case we have the dis- 
crete energy levels as schematically illustrated in Fig.[2j 
To calculate the escape rate T in Eq. (|5j) we employ the 
mean field approximation for x/j, that is, ip 2 and ip in 
Eq. ([6c]) are approximated with their expectation values. 
Then, at zero temperature we find (^)th = and 
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The finite value of (iJj 2 } originates from the zero-point 
motion of the "quantized" JL mode. Under this approx- 
imation we find the effective Lagrangian of single degree 



of freedom as follows, 
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where V^ m?e ff is the renormalized potential 
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Then, in this approximation the expectation value 
if({0},{0};/3) in Eq.© is reduced to K({0},{0};/3 = 

oo) - C)=o°° )= ° ™ exp(-ft- f«>=~ L^drl 
which can be evaluated in the standard instanton 
approximation—. Hence, the MQT escape rate corrected 
by the zero-point motion of ip is 
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where uj p (I) = uj p [(1 - e) 2 - J 2 ] 1 / 4 , V = 
h 2 [u P (I)} 2 cot 2 6>o/3^c, sin(9 = J, and I = I ex /I c . 

Figure [3] shows a contour map of the ratio T/Tq in the 
(J ex /J c vs. ojp/ojjl) -plane with To being the escape rate 
without correction, i.e., e = 0. It is seen that the escape 
rate is drastically enhanced in a wide parameter region. 
In particular, the enhancement is pronounced in the re- 
gion of large ojp/ojjl- As seen in Eqs. (j7a|) and (|7b|) . the 
Josephson coupling energy is renormalized by the zero 
point motion of ip and the renormalized one is decreased 
from the bare one since e > 0. As a result, the tunneling 
barrier for is lowered as schematically shown in Fig.[4j 
which causes the strong enhancement of the escape rate. 
In fact, R(e) = Vo/hup(I) is smaller than R(e = 0) for 
fixed I when < e < 1, that is, the exponent in Eq. flH]) 
is decreased. Thus, the renormalization increases T. It 
should be also noted that the zero-point fluctuation be- 
comes larger as the frequency of the JL mode decreases. 
Thus, the considerable enhancement of F occurs for the 
system with a lower value of cjjl- The MQT in the con- 
ventional systems is subject to the ratio Ej/Eq, which is 
an important parameter for designing a superconducting 
Josephson qubit 5 . In the system with multiple tunnel- 
ing channels the ratio co>p/cjjl(oc Ej / E- in ) also affects the 
characteristics of the MQT. 

In this paper we have focused on the tunneling process, 
\6 = 0, ip = 0) |0, 0) and clarified the effect of the JL 
mode on the MQT. We mention that such a process is 
not the unique one that contributes to the MQT rate in 
this system, since a system with two degrees of freedom 
generally have many tunneling routes. For example, the 
tunneling process in which the quantum switching in the 
6^ channel takes place successively after the switching 
in the 0^ channel will be also possible in the present sys- 
tem—. In this case the escape rate can be calculated from 



the transition process l^ 1 ) =0) — » |0) with 
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where f(t) is a time-dependent c- number function. This 
tunneling process is analogous to the MQT under a peri- 
odically time-dependent perturbation 25 . It is also noted 
that the relative phase difference ip might play a role of an 
environmental variable for through the term linear in ijj 
in Eq. (|6c|) . The MQT rate in this process can be evalu- 
ated, using the influential functional integral method—^. 
The competition between the zero-point fluctuation and 
the "dissipation" occurs in this case. The enhancement 
via the JL mode may be superior to the reduction from 
such dissipation when g+ > \g~\. 

We also mention that our theory for the MQT in the 
hetero Josephson junctions can be extended to the case of 
intrinsic Josephson junctions (IJJ's) with multiple tun- 
neling channels—. The MQT in such systems will be 
observed in several highly-anisotropic layered iron-based 
superconductors recently discovered— - —. In the IJJ's 
correction due to the JL mode for the corporative MQT 
among the junctions 3Q i 31 will be expected. 



Finally, we remark that the present theoretical predic- 
tion relies on the coexistence of the Josephson-plasma 
and the JL modes. Since observation of Leggett's mode 
in a bulk MgB 2 sample has been reported 32 and in junc- 
tions with MgB 2 33 i 34 , we expect that such a collective 
mode can be detected in a junction system and the the- 
ory will be verified experimentally. 

In summary, we have constructed a theory of the MQT 
in hetero Josephson junctions with multiple tunneling 
channels. We have clarified that the zero-point fluctu- 
ation of the relative phase differences brings about the 
drastic enhancement of the MQT escape rate. The en- 
hancement is gigantic when the JL mode has a lighter 
mass than that of the Josephson-plasma. 
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Science. 



1 A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 211 
(1981). 

2 R. Rajaraman, Solitons and Instantons: An Introduction 
to Solitons and Instantons in Quantum Field Theory (El- 
sevier, Amsterdam, 1987). 

3 E. Simanek, Inhomogeneous Superconductors: Granular 
and Quantum Effects (Oxford University Press, New York, 
1994) Chap.4. 

4 A. Larkin and A. Varkamov, Theory of Fluctuations in 
Superconductors Revised Edition (Oxford University Press, 
New York, 2005) Chap. 13. 

5 J. Q. You and F. Nori, Phys. Today 58, 42 (2005). 

6 J. Clarke and F. Wilhelm, Nature (London) 453, 1031 
(2008). 

7 M. Nakahara and T. Ohmi, Quantum Computing: From 
Linear Algebra To Physical Realizations (CRC Press, 
Talyor & Francis Group, Boca Raton, 2008) Chap. 15. 

8 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. 130, 3296 (2008). 

9 M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 
101, 107006 (2008). 

10 J. Paglione and R. L. Greene, Nat. Phys. 6, 645 (2010). 

11 J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, 
and J. Akimitsu, Nature (London) 410, 63 (2001). 

12 A. Brinkman and J. M. Rowell, Physica C 456, 188 (2007). 

13 X. X. Xi, Rep. Prog. Phys. 71, 116501 (2008). 

14 Y. Ota, M. Machida, T. Koyama, and H. Matsumoto, 
Phys. Rev. Lett. 102, 237003 (2009). 

15 T. Koyama, Y. Ota, M. Machida, Physica C 470, 1481 
(2010). 

16 Y. Ota, M. Machida, and T. Koyama, Phys. Rev. B 82, 
140509(R) (2010). 

17 D. F. Agterberg, E. Demler, and B. Janko, Phys. Rev. B 
66, 214507 (2002). 

18 A. J. Leggett, Prog. Theor. Phys. 36, 901 (1966). 

19 S. G. Sharapov, V. P. Gusynin, and H. Beck, Eur. Phys. 
J. B 30, 45 (2002). 

20 Y. Ota, M. Machida, T. Ko yama, and H. Aok i, Phys. Rev. 
B (to be published); e-print arXiv: 1008.32121 



21 H. Shimakage, K. Tsujimoto, Z. Wang, and M. Tonouchi, 
Supercond. Sci. Technol. 17, 1376 (2004). 

22 X. Zhang, Y. S. Oh, Y. Liu, L. Yan, K. H. Kim, R. L. 
Greene, and I. Takeuchi, Phys. Rev. Lett. 102, 147002 

(2009) . 

23 S. Schmidt, S. Doring, F. Schmidl, V. Grosse, P. Seidel, 
K. Iida, F. Kurth, S. Haindl, I. Monch, and B. Holzapfel, 
Appl. Phys. Lett. 97, 172504 (2010). 

24 Y. Chizaki, H. Kashiwaya, S. Kashiwaya, T. Koyama, and 
S. Kawabata (unpublished). 

25 M. P. A. Fisher, Phys. Rev. B 37, 75 (1988). 

26 S. Kawabata, T. Bauch, and T. Kato, Phys. Rev. B 80 
174513 (2009). 

27 H. Ogino, Y. Matsumura, Y. Katsura, K. Ushiyama, S. 
Horii, K. Kishio, and J. Shimoyama, Supercond. Sci. Tech- 
nol. 22 075008 (2009). 

28 H. Nakamura, M. Machida, T. Koyama, and N. Hamada, 
J. Phys. Soc. Jpn. 78, 123712 (2008). 

29 H. Kashiwaya, K. Shirai, T. Matsumoto, H. Shibata, H. 
Kambara, M. Ishikado, H. Eisaki, A. Iyo, S. Shamoto, I. 
Kurosawa, and S. Kashiwaya, Appl. Phys. Lett. 96, 202504 

(2010) . 

30 M. Machida and T. Koyama, Supercond. Sci. Technol. 20, 
S23 (2007). 

31 S. Savel'ev, A. O. Sboychakov, A. L. Rakhmanov, and F. 
Nori, Phys. Rev. B 77, 014509 (2008). 

32 G. Blumberg, A. Mialitsin, B. S. Dennis, M. V. Klein, N. 
D. Zhigadlo, and J. Karpinski, Phys. Rev. Lett. 99, 227002 
(2007). 

33 Ya. G. Ponomarev, S. A. Kuzmicheva, M. G. Mikheeva, M. 
V. Sudakovaa, S. N. Tchesnokova, N. Z. Timergaleeva, A. 
V. Yarigina, E. G. Maksimovb, S. I. Krasnosvobodtsevb, 

A. V. Varlashkinb, M. A. Heine, G. Mullerc, H. Pielc, L. G. 
Se vasty anovad, O. V. Kravchenkod, K. P. Burdinad, and 

B. M. Bulychevd, Solid State Commun. 129, 85 (2004). 

34 A. Brinkman, S. H. W. van der Ploeg, A. A. Golubov, H. 
Rogalla, T. H. Kim, and J. S. Moodera, J. Phys. Chem. 
Soilds 67, 407 (2006). 



